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Abstract 

Quantum fluctuations for a massless scalar field in the background metric of spherical 
impulsive gravitational waves through Minkowski and de Sitter spaces are investigated. It 
is shown that there exist finite fluctuations for de Sitter space. 
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Introduction 



The possibility that at the root of the mechanism explaining galaxy formation cosmic 
strings lie is not still ruled out I 1 . We need further experimental data on the anisotropies 
in the cosmic microwave radiation, lensing of quasar images, or gravitational radiation 
stemming from the decay of strings to accept or reject this alternative to inflationary 
quantum fluctuations Z 1 ' 2 ' 3 . 

However, besides their cosmological implications cosmic strings are also of interest 
from the point of quantum field theory in curved space-time. Since the space-time gener- 
ated by a straight cosmic string possesses a conical structure, it gives rise to finite vacuum 
fluctuations of quantized fields. These effects have been calculated by many authors 4 ' 5 . 
Another effect is seen when the straight cosmic string splits into two ends. This results in 
the emission of a spherical wave. 

Exact solutions of Einstein's field equations for spherical impulsive waves emitted by 
snapping cosmic strings was found by Nutku and Penrose / 6 . It is given by the metric 

dsi=(2dudv + 2u*\d( + vG{v)hm f 
\ 2u 

where h is the Svhwarzian derivative of an arbitrary function /(C)- C is t ne stereographic 
coordinate on the sphere, u is a Bondi-type luminosity distance and v is a null coordinate 
which can be regarded as retarded time. 

One may investigate whether the spherical impulsive wave emitted when a cosmic 
string snaps gives rise to finite vacuum fluctuation phenomena. In Ref. 7 we investigated 
the vacuum fluctuations in the background metric of special cases of these solutions. 

One can construct similar solutions for impulsive gravitational waves propagating 
through the de Sitter universe / 8 . Here one has to simply multiply the Minkowski solution 
by a conformal factor (l + ^p)~ 2 . Although the Nutku- Penrose solution is not conformally 
related to Minkowski space., an obstruction to the presence of finite vacuum fluctuations 
may sneak in since we are perturbing about the Minkowski solution. This new solution in 
de Sitter space has a parameter, the cosmological constant A, with the dimensions of mass 
squared. In thuis case, perturbation theory is performed about a space where conformal 
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symmetry is explicitly broken. This allows us to investigate if our null results in Ref. 7 
persists in the absence of conformal symmetry even in the space about which we perform 
our perturbation expansion. 

To search for vacuum fluctuations, the usual method is the calculation of vacuum 
expectation value, VEV, of the the stress-energy tensor of a scalar field in the background 
metric proposed for the spherical wave / 9 . This is achieved by calculating the two point 

5(x x') 

function, Gf(x, x') for the model as a solution of the equation (D + £R)Gf = — j-J an< ^ 
differentiating the result to obtain < 0|T Mi ,|0 > after the coincidence limit is taken. The 
d'Alembertian, denoted by □, in the above expression is written in the background metric 
of Eq.(l). R is the Ricci scalar and £ is a constant which determines the coupling of the 
field to the metric. 

Here we will first perform the calculation by expanding about the Minkowski solution. 
This means in calculating VEV we will use the vacuum for the empty Minkowski space- 
time with no boundaries. In other words we will define our vacuum state in the usual way 
by choosing positive frequency modes with respect the Killing vector of flat Minkowski 
space-time. Then we will obtain the de Sitter solutions by multiplying this solution by 
conformal factors. We will obtain the de Sitter solution by multiplying this solution by 
conformal factors . 

For technical reasons we first take a massive scalar field, but take the mass equal to 
zero at the end of the calculation not to introduce any mass parameters to the Minkowski 
calculation. We have more control on the consequences of such a parameter when it is 
introduced via the de Sitter universe construction. We use the conformal coupling of the 
scalar field to the metric. This is important in the de Sitter case where the Ricci scalar is 
not zero. In the Minkowski case the Ricci scalar is zero; so, whether we use minimal or 
conformal coupling does not matter. 

In Ref. 7 we studied a special case of the metric given in Eq. (1) and used the form 
/ = (C) 1+S+ie - The arbitrary function in the solution will be the Schwarzian derivative if 
this function. 

We took only the first order terms and found that h is proportional to 5 and e. Since 
in the Nutku-Penrose solution the impulsive wave is generated in a process where a cosmic 
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string snaps, we want a small parameter multiplying any possible non zero fluctuation 
effect to act as a signature of the string. Here 5 and e serve this purpose. 

In this paper we take a generalized form of this function as / = J 
took this form for / , since we knoe that the Schwarzian derivative of this expression is 
zero when e is zero, h is proportional to e in the first order. 

Our previous choice for / is a degenerate form of this expression, where either A and 
D or C and B are taken to be zero. We take 5 zero in the new calculation since we know 
from our previous work I 1 that for fluctuations the case when 5 is not equal to zero looks 
essentially the same as the case when e is not to zero. The singular part of the Greens 
function is the same, only factors multiplying this term changes; so, we conclude that 
taking only one of the parameters finite suffices. 

If we agree not to take both A and C in the above expression equal to zero, we can 
reduce the new / to the form / = I p ) up to an overall factor in front of this 



term which can be absorbed in the expansion parameter e. In our specific calculation we 
take ^ = — 1, jj = 1 for convenience. Since we studied the A = and the C = cases 
seperately, we think that our new case exhausts all the interesting cases. Whether there 
are vacuum fluctuations or not should not depend upon the location of the zeroes and poles 
of /; so, taking the poles at plus and minus one should be able to represent the general 
behaviour of this set of trial functions. 

In the main part of the article, we describe the stress-energy tensor calculation when 
/ = ( ^pj J • In Section III we apply our results to the de Sitter case. Since the de 

Sitter solution is just conformally related to the Minkowski case, we obtain the de Sitter 
space Greens functions from the Minkowski ones just by multiplying them by conformal 
factors. The end results is different, though. In all the cases studied, we find one non- 
vanishing component of the vacuum expectation value of T vv which is proportional to 
8(v) A 2 . The metric we used had only one non identically vanishing component for the 
curvature tensor which is also proportional to the Dirac delta function of the variable 
v .The same behaviour is seen in the stress-energy tensor in de Sitter universe. 
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Calculationg for the Wave Propagating in Minkowski Space 

We start with the Nutku- Penrose metric, given in eq.(l) for the case where h is the 

//_l\ 1+iv/ 2 £ 

Schwarzian derivative of / = ( j and treat only first order terms in e. Here taking 

complex values for the exponent corresponds to taking a rotating string, as explained in 
Reference 6. 

In the first order in e,h is given by 

h ^ C (2) 

~ (x + iy + l) 2 (x + iy - l) 2 

f° r C = (2)1/2 • This choice gives us the metric 

ds 2 = 2dudv - (u 2 + eP^dx 2 - (u 2 - eP^dy 2 + 2eP 2 dxdy (3) 

where 

( X 2 — y 2 — I) 

Pi = 64uvxy± ^ " ( 4 ) 

[(x 2 - y 2 - l) 2 - Ax 2 y 2 ] 
P 2 = I6uv ^ ( 5 ) 

P= ( X 2 -y 2 -if + 4 X 2 y 2 . (6) 

Just note that for technical purposes we first do the calculation for the massive case. 
Here the mass used can be taken as an infrared parameter. We take the zero mass limit 
before the coincidence limit is taken. We derive the d'Alembertian operator (— ^) _1 / 2 L, 
written for the Nutku-Penrose metric. The operator, L , reads, for v > 0, 



d 2 d d 2 d 2 

2u 2 —— + 2u— - -— — — — ?r -\- 2m 2 u 2 + 



d 2 d 2 , „ d 2 1 



^ife-7or)-2P 5 



dx 2 dy 2 ' 2 dxdy 



(7) 



dudv dv dx 2 dy 2 u 2 
We note that the Nutku-Penrose metric does not allow first order derivatives in x and y, 
for any h used in the original metric in this order. 

In order to solve Eq.(7) we use a perturbative scheme. The zeroth order part, which 
we call Lq, is the flat Minkowski space d'Alembertian function written in our coordinate 
system. Our vacuum will be that of the Minkowski space with no boundaries. 

At this point we may use an expansion of the Feynman Green Function using the 
integral equation approach. Then 

Gf = G F + G F VG F + 
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where G F is the free Feynman propagator and V is the perturbing part. 

We will, instead, start with the Sturm-Liouville equation L(p\ = and the explicit 
construction of the Green function Gf will use the eigenfunctions of this equation. Gf 
will be given by 

<j>x{x)P x {x') 



Gf =-y: 



A 



(8) 



Here the eigenfunctions of the Sturm-Liouville problem form a complete set. We perform 
the calculation to first order in e. 

We expand the operator L, the eigenfunction <p\ and the eigenvalue A, in powers of e, 



which gives 



(L + eLi)(0 o + e0i + • • •) = (A + eAi + • • -)(0o + efa + • • •)• 

L (fio = A o 0o, 
Li0 o + ^o^i = Ai0 o + A o 0i 



where 



L = (2u 2 



2 



dudv 



+ 2u— 



d d 2 



d 2 



2„.2 



dv dx 2 dy 2 



) + 2m z u 



d 2 



r)-2P 2 - 



(9) 

(10) 
(11) 

(12) 
(13) 



<9x 2 <9|/ 2 7 ^ dxdy 
L is the d'Alembertian for the Minkowski case. Pi and P 2 are as defined in Eq.s (4) and 
(5). 

The zeroth order solutions 0o and Aq are found easily . 
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(14) 



{2n) 2 u^2\R\ 

Here A = K — k\ — k\. ki, k 2 , R and K are constants definig the different modes. Note 
that 0o is also a solution in Minkowski space. 

We also find Ai = (0o, £i0o) = due to the particular form of the operator L\. 

To solve the inhomogenous equation for 0i, we take <f>i = 4>of. Then our eq.(ll) 
reduces into L±f = vH, explicitly written as 



rt .x, 2 d n . n d d , d 2 d 2 
2iRu 2 — - 2i(k 1 — + k 2 —) - — - — 
du ox ay dx z dy z 



+ 2u' 



d 2 



+ 



%K d 



dudv R dv 



f 
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V 

u 



64xy(x 2 -y 2 - l){k\ - k 2 2 ) hk^x 2 - y 2 - l) 2 - Ax 2 y 2 ] 



[(x 2 - y 2 - l) 2 + 4x 2 y 2 ] 2 ' [(x 2 - y 2 - l) 2 + 4x 2 y 2 ] 2 
A particular solution of this equation is of the form 



(15) 



/ = vF^x.y.u) + F 2 {x,y,u). 



(16) 



We find this ansatz by assuming the solution to be of the form 

/ = Gi(w)Fi(x, y, u) + F 2 (x, y, u). 

When this form is substituted in Eq. (15), the form given in Eq.(16) emerges, with two 
new coupled differential equations for F\ and F 2 . 



L 2 F\ = H, 



L 2 F 2 = (-2u 2 - 



(17) 



(18) 



du R 

Here H is as defined inside the square parenthesis on the righr hand side of eq. (15) and 
L 2 is given as 

(19) 



n . n 2 9 n . n d , d . d 2 d 2 
L 2 = 2iRu 2 — - 2i(k 1 — + k 2 — ) -jr^-jr^ 
ou ox ay ox z ay 2 



Since 0i is the product of <po with /, the boundary conditions are dictated by <fi$, and 
the form given by Eq.(16) obeys the boundary conditions fixed by the zeroth order term, 
which is the solution in empty Minkowski space. 

At this point we change the variables and define s = ^. We use both variables in our 
expressions. We also note that if we use the new variables z = x + iy,z = x — iy, both the 
inhomogenous term H and the operator itself seperate. We find that we get 



d ( d d \ d 2 

-2 iR --2iUk 1 + ik 2 )- + (k 1 -ik 2 )-)-,— 



Fi = 



-8is 



(ki-ik 2 y 



(fci + ik 2 ) : 



(z-l) 2 (z + l) 2 (z-l) 2 (z + l) 2 _ 
and a similar equation for F 2 . 

Upon solving the new differential equations, we find particular solutions 



(20) 



and 



where 



F x =s[2iA;i(tan- 1 2y „ - (- - + -)) 

1 K x 2 + y 2 -l y (x-l) 2 + y 2 (x + l) 2 +y 2/; 

_ , (x + l) 2 + y 2 2(x + l) 2(x-l) „ 

2 ^° S (x-l) 2 + y 2 (x + l) 2 + y 2 (x-l) 2 +y 2 ^ ( ' 

t,/ -i % y, (^-i) 2 + ?/ 2 n 



" 2= ^ ~2R ~ 7 Ml(x '^ + 2(fc 2 + k 2 ) ^ klM ^ x ^ ~ k 2M 3 (x,y)), (22) 



Ml (x,y) = -2ylog ( ^^tya ) " 2ajtan_1 ^1 + 2Xtan_1 ^1 (23) 
M 2 (x,y) = 4*, log +2(x 2 -, 2 -l)(tan- 1 ^"tan" 1 ^)+% (24) 

M 3 (x, y) = (2* 2 -2y 2 -2) log ffi "^(tan" 1 ^y-tan" 1 ^0+4* (25) 

Now these solutions are used to form the Greens function, Gf, of the original operator 
given in eq.(7). Since the zeroth order part of is the same as the flat space Greens 
function, we calculate only the first order part. 
To obtain Gf, in first order, we take 

* r r *. r iK r at^^^**^,^-*) 



G F = - dh dk 2 dK tt^^^- a e l 

2 V-oo J -00 J -00 \R\ (2tt)W 

K-fc 2 -fc 2 • (2bj 

with f(x,y,u,v) defined in Eq. (16). Here we use the Schwinger prescription to regular- 
ize the denominator and use the formula A ^_ iS = i J °° e~ iaA ~ 5a da where 8 is a positive 
constant approaching zero. The infinite integration over ki,k 2 is regularized by the de- 
nominator. The result is the Feynman propagator with its usual contour. 

We need to calculate the vacuum expectation value of the stress-energy tensor 
T M j,. We are particularly interested in whether < T vv > has a finite part. Since in T vv 
x or y derivatives do not exist, we take the coincidence limit in x and y at this stage. 



This puts terms of the type (x — x') 2 and (y — y') 2 to zero. Such terms either multiply 
the expressions obtained, or appear in the 'geodesic distance' term a = (u — u'){v — v') — 
^-((x — x') 2 + (y — y') 2 ), which appears in the denominator. Taking the coincidence limit 
in x and y only simplifies the calculation at this stage and does not affect the end result. 

When we take x equals x' and y equals y' ', all the terms that are odd in k\ and k 2 
vanish. The remaining terms are given by 

GF(s,s'v,v'x,y) = j— — — r / [ dK [ dk\ [ dk 2 [ dae~^ nr ^~~^ 

1 ' y> uu'2{2^YJ 2\R\ J J J Jo 

x e iR(v-v') e im 2 e -ia (K-k*-k>) -aS 

x [ ^( S e( v ) - s 'Q{v')) + i{Q{v) + e(v')fj M 1 {x, y) + 2i?M 4 (x, y)(vQ(v) - v'Q{v'))] 

(28) 

where 

/ \ 1 , (x + 1) 2 + y 2 i 2y , nn . 

M 4 (x. y) = -y log L _ > y - - .tan-' ^^-p (29) 

and M 1 (x,y) was defined in Eq. (23). 

We perform the integrations over k\,k 2 ,K and R. The result of the a integration 
can be shown to result in Hankel functions which degenerates into a monomial when the 
m going to zero limit is taken. Hence for the massless field case the Greens Function Gf 
reads 

V16tt/ l (w-m')(^-^ ) V (u-u') K KJ y ") V ' 

\^_ t;/ ) - )] (3°) 

We note that the singularity structure in the coincidence limit of this expression 
displays exactly that of the free Minkowski Green Function. There is just a modulating 
factor which is finite in the coincidence limit. Since our perturbative solution (pi is given as 
a factor times the Minkowski solution, the singularity structure of the Minkowski vacuum 
is carried over this solution. This part may be regularized by renormalizing the free 
parameters of the theory. 
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The empty space Green function goes as Here a is the square of the geodesic 

distance between the two points given by Eq. (27). A well known theorem / 10 states that 
if the metric is of the Minkowskian shape for a certain region of space-time, the Green 



where A,B,C are finite quantities at the coincidence limit; a is the geodesic distance be- 
tween the points of Gf- Our hope, in this calculation, was to find a finite term C given 
in the above formula. Unfortunately our end had B and C identically zero. We choose 
our boundary conditions of the renormalized Green function of the renormalized Greens 
function so that we get zero as u and v go to infinity. This condition fixes Cq in the 
expansion C = ECjctj, therefore completely fix C uniquely equal to zero. 

Another way of obtaining this result will be applying the Adler-Lieberman-Ng pre- 
scription / 12 to this problem. An ambiguity may arise only by a local conserved tensor 
which comes out to be proportional to a<i of the adiabatic expansion in the DeWitt- 
Schwinger method Z 13 . This term is identically zero in our case; so, even this possible 
ambiguity is not present. Other references to the uniqueness property of the VEV of the 
stress-energy tensor can be found in Ref. 14. 

Calculation for the Wave Propagating in de Sitter Universe 

We conjecture that we have to introduce a mass parameter to our space in order to 
get non zero result for the vacuum fluctuations if we insist performing our calculation 
perturbatively. One can check whether this is true by doing the similar calculation in de 
Sitter universe. Indeed we find that our end result changes if we look at the spherical 
impulsive gravitational waves propagating through the de Sitter universe. It is shown that 
/ 8 an impulsive spherical gravitational wave solution exists in de Sitter space. Here one 
has to simply multiply the Minkowski solution by a confer mal factor (1 + ; ^p) where A 
is the cosmological constant which comes out proportional to Ricci scalar in this metric. 
The metric now reads 



Function will be of the Hadamard type ' n , 



G F = — + Bloga + C, 



a 




(31) 
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We take the same h as before, as given in eq.(2). Since we use conformal coupling, the 
operator L is modified to include also the Ricci scalar term for this metric which is no 
longer zero. 

We know from general arguments that ^ 9 

G| = (l+^)Gf(x,x')(l + ^) (32) 

where Gp and Gp are the de Sitter and Minkowski space Greens functions respectively. 
Note that 

(i + + ^) = (i + + A ( „ _ K > - </) 

+(A)» (-(„ - u>?V> - („ - v>?V> + (U -"' ) >-' /)2 ) (33) 

where [/ = U = Using this expansion we see that for Gp we get terms in the 

coincidence limit that give finite terms to < T vv > 



lim 



d 



2 



u—>u',v—>v' dvdv' 



A 2 (u-u') 2 (v-v r 



2 



G%{x,x') 



(24)2 

= wf, ( 2 — 1 iTi - 2 * tan_1 " %log (IwtS) (34) 

We find that in first order calculation for spherical impulsive waves propagating 
through the de Sitter universe, we get one non-zero component of the stress-energy tensor. 

If we perform the similar calculation for the cases studied in Ref.7, we get nonzero 
contribution, given in the same way. Here h in eq.(l) comes from / = (Q 1+5+ie - Using 
the results given in Ref.7, we get 

-A 2 u5(v) 



T = 

mi 



(5\og{x 2 + y 2 ) -4etan" 1 |) . (35) 



2(96)% 

This result is finite, contrary to the case for propagating in Minkowski space. 
Conclusion 

Here we calculated the stress-energy tensor of a scalar field in the space-time of a 
spherical impulsive gravitational wave propagating through the Minkowski and de Sitter 
universes. 
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We verified the common knowledge that fluctuations which are null in Minkowski 
space may become finite in de Sitter space. We found that VEV of one component of the 
energy momentum tensor in first order is proportional to 8(v) which is the signature of 
the impulsive wave solution, in de Sitter space. This is the same factor to which the only 
nonzero component of the curvature tensor is also proportional. 

Our results are found in first order perturbation theory, but are proportional to the 
square of the curvature scalar. They are an example of extracting non trivial information 
out of first order perturbation theory. 

As a technical remark note that first order solutions to the Sturm Liouville problem 
arising from this metric can be expressed in terms of sum of holomorphic and antiholo- 
morphic functions. This occurs after we factorize the solution into two parts, one part 
proportional to the flat space solution. This fact simplifies our calculation, since then sec- 
ond order differential equation decomposes in to a pair of first order equations which are 
easily integrated. In our previous papers / 7 , we used brute force Greens function method 
to integrate the inhomogenous equation. Since we have to perform integrals where the 
integration range is infinite, there were points where regularizations were used. In the new 
method we do not need any regularizations in calculating Gp- The result is free from any 
possible ambigiuties the regularization may cause. 

We reduce first order eigenfunction calculation to the solution of a couple of first 
order differential equations. This is possible since both the inhomogenous term and the 
differential operator decomposes into a sum of two terms if a set of variables are used. 
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